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We study a particle production at the collision of two domain walls in 5-diniensional Minkowski 
spacetime. This may provide the reheating mechanism of an ekpyrotic (or cyclic) brane universe, 
in which two BPS branes collide and evolve into a hot big bang universe. We evaluate a production 
rate of particles confined to the domain wall. The energy density of created particles is given as 
p Ri 20g Nh m^ where g is a coupling constant of particles to a domain-wall scalar field, A^f, is the 
number of bounces at the collision and m,, is a fundamental mass scale of the domain wall. It does 
not depend on the width d of the domain wall, although the typical energy scale of created particles 
is given by tj ~ 1/d. The reheating temperature is evaluated as Tr ~ 0.88 g Nf^ . In order to 
have the baryogenesis at the electro-weak energy scale, the fundamental mass scale is constrained 
as m^>l.l X 10^ GeV for g - 10"^ 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

The big bang theory is very successful because it natu- 
rally explains the evolution of our universe from the nu- 
cleosynthesis to the present time with many observational 
data. However, it contains some theoretical key problems 
such as the flatness and the horizon problems[l|, Ql- So 
far, only the idea of inflation provides a resolution of 
those problems. Not only it gives some picture of the 
earlier stage of the universe before a big bang but also 
it seems to be supported by some of recent observational 
data on CMB. While, it is still unclear what is the ori- 
gin of inflaton. So far, there is no convincing link with 
fundamental unifled theories such as a string/M theory. 

Recently a new paradigm on the early universe has 
been proposed, which is so-called a brane world jH 0|- 
Such speculation has been inspired by recent develop- 
ments in string/M-theory[3, M Q- There has been 
tremendous works in this scheme of dimensional reduc- 
tion where ordinary matter flelds are confined to a lower- 
dimensional hypersurface, while only gravitational flelds 
propagate throughout all of spacetime. In particular, 
it was shown that the 10-dimensional Es x Eg het- 
erotic string theory, which is a strong candidate to de- 
scribe our real world, is equivalent to an 11-dimensional 
M theory compactifled to M^° x SV-^2 0- Then 
the 10-dimensional spacetime is expected to be com- 
pactifled into M^ X CY^ where M" and CY^ are 
4-dimensional Minkowski spacetime and 6-dimcnsional 
Calabi-Yau space, respectively. Randall and Sundrum 8] 
also proposed a new model where four-dimensional New- 
tonian gravity is recovered at low energies even without 
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compact extra dimensions. Based on such a new world 
picture, inany cosmological scenario have been studied 
E E Il3- See also recent reviews 51 El El El. We 
have found some deviations from standard cosmology 
by modiflcations of 4-dimensional Einstein equations on 
the brane [ifij. even for the case there is a scalar fleld in 
bulk[l3. 

In such a brane world scenario, for resolving the above- 
mentioned theoretical key problems in the big bang the- 
ory, new idea of the early universe has been proposed, 
which is called the ekpyrotic scenario or a cyclic uni- 
verse scenario [l^ E^ ■ It is based on a collision of two 
cold branes. The universe starts with a cold, empty and 
nearly BPS ground state of M theory, which contains 
two parallel branes at rest. Two branes approach each 
other and then collide. The energy is dissipated on the 
brane and the big bang universe will start. The BPS 
state is required in order to remain a supersymmetry in 
a low-energy 4-dimensional effective action. The visi- 
ble and hidden branes are flat, which are described by 
a Minkowski spacetime, but the bulk is warped along 
the flfth dimension. Since this scenario is not only moti- 
vated by the fundamental unifled theory but also it may 
resolve the theoretical key problems such as the flatness 
and horizon problems, it would be very attractive. There 
are also many discussions about density perturbations to 
see whether this scenario is really a reliable model for the 
early universelH IH El HI . 

On the other hand, even though there are some works 
by 24J, the reheating process itself in this scenario has 
not been so far investigated in detail. Hence, in this 
paper, we study how we can recover the hot big bang 
universe after the collision of the branes. Here we inves- 
tigate quantum creation of particles, which are conflned 
to the brane, at the collision of two branes. It may be 
difflcult to deal properly with the collision of two branes 
in a basic string theory. Hence, in this paper, we adopt a 



domain wall which is constructed by some scalar field as 
a brane, and analyze the collision of two domain walls in 
a 5-dimensional bulk spacetime. Some works also adopt 
such a picture [23, llg ■ It is worth noting that there is a 
thick domain wall model for a brane world|27J. 

In order to analyze particle creation at the brane colli- 
sion, in this paper we consider the simplest situation. We 
discuss that two domain walls collide in a 5D Minkowski 
spacetime. In Sec. m we analyze the collision of two 
domain walls. Then, in Sec. IIIII we investigate particle 
creation on the wall at the collision. Applying the par- 
ticle production to the energy dissipation of the brane, 
we discuss the reheating mechanism of a brane universe. 
We use the unit oi c — h= 1. 



II. COLLISION OF TWO DOMAIN WALLS 

A. Basic equations and initial setting 

We study the collision of two flat domain walls in 5- 
dimensional (5D) Minkowski spacetime. To construct a 
domain wall structure, we adopt a real scalar field $ with 
a double-well potential. 
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where the potential minima are located at <I> = ±77. 

We discuss the collision of two parallel domain walls. 
The scalar field is assumed to depend only on a time 
coordinate t and one spatial coordinate y. The rest three 
spatial coordinates are denoted by x. We rescale the 
parameters and variables by rj (or its mass scale rriri — 
rf/^) as 



t — m„t , y — rriny ,<&=—, A = rrinX 



(2.2) 



In what follows, we omit the tilde in dimensionless vari- 
ables for brevity. 

The equation of motion for $ in 5D is given by 



$-$" + A<I>($2- 1) = 0, 



(2.3) 



where ' and ' denote d/dt and d/dy, respectively. 

Eq. H2.3(l has not only two stable vacuum solutions 
^± — ±1 but also a static kink solution (K), which is 
topologically stable. The latter one is called a domain 
wall, which is described by 



$K(y) = tanh 



(2.4) 



where d — a/2/A is the thickness of the wall '2q . The 
antikink solution {K) is obtained from Eq. (|2.4() by re- 
flecting the spatial coordinate y as ^niy) = ^K{—y) = 
—^K{y)- A single domain wall solution moving with a 



constant speed v in the y direction is given by boost- 
ing the solution (|2.4|l . We find a moving domain wall 
solution, whose initial position is located at y = 0, as 



^v{y,t) = tanh 
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(2.5) 



where 7 = 1/Vl — w^ is the Lorentz factor. 

In order to discuss the collision of two domain walls, we 
first have to set up the initial data. We put a kink and an 
antikink solution, which are separated at large distance, 
and are moving toward each other with velocities v and 
—V. The explicit initial configuration is given by 

$(2/, 0) = $,(2/ + j/o, 0) - $_,(y - 2/0, 0) - 1 , (2.6) 

where y = ±2/0 is the initial position of the walls. The 
spatial separation of two walls is given by 2j/o, and as long 
as the separation distance is much larger than the thick- 
ness of a wall (j/o ^ d), the initial profile 1)2. 6|l gives a 
good approximate solution for two moving domain walls. 
The initial value of $ is also given by 



$(2/, 0) = $,(2; + 2/0, 0) - ^-^,iy ~ 2/0, 0) 



(2.7) 



Using this initial values, we solve the dynamical equation 
(|2.3() numerically, whose results will be shown in the next 
subsection. 



B. Time Evolution of Domain Walls 

We use a numerical approach to solve the equations 
for the colliding domain walls. The numerical method is 
shown in Appendix ^ We have two free parameters in 
our simulation of the two-wall collision, i.e. a thickness 
of the wall d = a/2/A and an initial velocity of the walls 
V. The collision of two walls was discussed by |23| in 4- 
dimensional Minkowski space. Although we discuss the 
domain wall collision in 5-dimensional Minkowski space, 
our basic equations are exactly the same as their case, 
and then we find the same results as theirs. In particular, 
the results are very sensitive to the initial velocity v. 

First we show the numerical results for two typical 
initial velocities, i.e. w = 0.2 and 0.4 in Figs. ^[21 The 
evolution of the energy density of $ is shown in Figs. Q] 
and 121 The energy density of the scalar field is given by 
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From Figs. ^ and El we find some peaks of the energy 
density, by which we can define the positions of moving 
walls (2/ = ±yw{t)). 

From Figs. ^andEl we see the detail of the collision 
as follows. For the case of initial velocity v = 0.4, 
we find that, the collision occurs once, while it does 
twice for the case of w = 0.2. To be precise, in the 
latter case, after two walls collide first, they bounce, re- 
cede to a finite distance, and then return to collide again. 




FIG. 1: Collision of two domain walls for the case of the initial 
velocity v — 0.4. The time evolution of the energy density of 
scalar field $ is shown from t = to 150. The maximum 
point of Pi, defines the position of a wall {y = j/w(i)). 




FIG. 2: Collision of two domain walls for the case of the initial 
velocity v = 0.2. The time evolution of the energy density of 
the scalar field <I> is shown from t = to 150. We find that 
collision occurs twice at f « 58 and 77. We set A = 1.0. From 
this figure, we find clearly that collision occurs twice. 



As shown by several authors [23, 120, IM; |33| , however, 
the result highly depends on the incident velocity v . In 
Appendix B, we show our detail analysis, which confirms 
the previous works. For a sufficiently large velocity, it 
is expected that a kink and antikink will just bounce off 
once because there is no time to exchange the energy 
during the coUisional process. In fact, it was shown in 
|29j that two walls just bounce off once for w>0.25. For 
smaller velocity, we find multiple bounces when they col- 
lide. The number of bounce during the collision depends 
complicatedly on the incident velocity. For example, the 
bounce occurs once for v = 0.4, while twice for v = 0.2. 
We also find many bounce solutions for other incident ve- 



locities as shown in Appendix B (see also |23). We find 
that the number of bounce is not monotonic function off, 
but is much more complicated. If we change the incident 
velocity slightly, the number of bounce changes. The ex- 
istence of a fractal structure in the parameter space is 
shown in Fig. 6 of |29l |. 



III. PARTICLE PRODUCTION ON A MOVING 
DOMAIN WALL 

A. Quantization of a particle on the domain wall 
and its production rate 

Once we find the solution of colliding domain walls, we 
can evaluate the time evolution of a scalar field on the 
domain wall. Since we assume that we are living on one 
domain wall, we are interested in production of a parti- 
cle confined to the domain wall. We assume that there 
is some coupling between a 5D scalar field $ which is 
responsible for the domain wall and a particle on the do- 
main wall. Because the value of the scalar field changes 
with time, we expect quantum particle production oc- 
curs. This may be important for a reheating mechanism 
of the colliding domain walls. 

Hence we have to know the value of the scalar filed <i> 
on the domain wall, i.e. $w(''') = ^{t,yw{t))- Since the 
wall is moving in a 5D Minkowski space, we have to use 
the proper time t of the wall, which is given by 

r = J dt^l - y^it) , (3.1) 

when we estimate the particle production in our 4- 
dimensional domain wall. 

We consider a particle on the domain wall which is de- 
scribed by a scalar field ij^. We assume that it is confined 
to the domain wall and couples to the scalar field $ as 



'■^'{r)^' 



(3.2) 



where g is a coupling constant. Here (J){t) is the 4- 
dimensional scalar field, which is induced from the 5- 
dimensional $ as 



(3.3) 



2d $w(-r) 



We have also assumed that the effective width of the do- 
main wall when it collides is about 2ci, which is confirmed 
by numerical simulation. In Fig. |21 we depict the spatial 
distribution of the scalar field $ when the domain walls 
collide. Note that <i>w(''") is the maximum amplitude in 
this distribution. 

The basic equation for ijj is then 
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f = 2dg^ = 



2%/2 
Ai72- 



(3.4) 



(3.5) 



as 



- 






- 


- 






- 


- 




1 y^ 

1 / 


- 


- 


\ 1 


1 / 


- 


- 


\l 


1/ 


- 


- 


\ 


1 


- 


- 


A 


i 


- 


- 


|\ 


\\ 


- 


- 


1 \ 


1 \ 


- 


- 






- 



FIG. 3: The spatial distribution of the scalar field <l> when 
the domain walls collide, d is the thickness of the wall. The 
scalar field has non-vanishing value for the effective width of 
2d at the collision. 



Since our background spacetime is the 4-diniensional 
Minkowski space, it is easy to quantize the scalar field 
■0. In a canonical quantization scheme |3.ll | . we expand V' 
as 

V'(T,a;) = ^[afcVfc(T)ufc(a;) + 4Vfc(T)w^(a;)] , (3.6) 



where Uk{x) — {2n) 2e'^'^'^. The wave equation H3.4|l 
for each mode is now 



v;fe + [fc'+r*wM]V'fe = o, 



(3.7) 



where k = \k\. Since two domain walls are initially far 
away each other, the value of $w is almost zero. We can 
quantize ip as a usual quantization scheme. The eigen 
function with a positive frequency is given by 



/ (in) 
V'fc 



1 



(3.8) 



where cjk — \/k'^ + g^^^iO) ~ k. We impose the equal 
time commutation relation for the operators a^ and aj, 

[afe,afe'] = 0, (3.9) 

(3.10) 
(3.11) 



[4,4,] -0, 



[ak,a\'\ = ^kk' , 



where aj; and a\. are an annihilation and a creation op- 
erators. We then define a vacuum |0 >i„ at r = by a^ 



afe|0>in=0 



V/c 



(3.12) 



After the collision of domain walls, we expect that the 
value of <i>vv again approaches zero (see the detail in the 
next subsection). We can also define the vacuum state 
|0 >out7 which is different from the initial vacuum state 
|0 >in. The eigen function of ^fc for r ^ oo is then given 



by a linear combination of i/jj, and ij)^ 
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(3.13) 

and the annihilation and creation operators as 

Ofe = akttk + Pla\ ' (3-14) 

where afe, Pk are the Bogolubov coefficients, which sat- 
isfy the normalization condition 



|afep-|/3fc|2 = l. 
The Hamiltonian of this system is given by 



(3.15) 
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(3.16) 



where : : is the normal ordering operation. The creation 
of the particles with mode k is evaluated as 



< 0|i„ : Hk:\0 >in= |/?/c|^Wfc as r ^ cx) 



(3.17) 



As a result, the number density and energy density of 
produced particles are given by 



n = j\l3k\^<fk, 

\h?ojkd'k 



(3.18) 
(3.19) 



B. Time evolution of a scalar field on the domain 
wall and particle production 

Now we estimate the particle production by the do- 
main wall collision. In Figs. 0]and[Sl we depict the time 
evolution of $w on one moving wall with respect to r. 
In Fig. n we find one collision point, which corresponds 
to a spike in Fig. ^ and two-bounce in Fig. |21 gives two 
spikes in Fig. |S1 

We also show the results for different values of the cou- 
pling constant A in Figs. |Sland[71(A = 10). If Fig. Elwe 
find that when A gets larger as A = 10, the spike of <I>w 
becomes sharp. The same thing happens in the case of 
two bounce (see Fig. [TJ . 

In Figs. I4I7I we find that <I>w begins to oscillate after 
the collision. We also find that the period of these os- 
cillations in Figs. El and [3 is shorter that those in Figs. 
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FIG. 4: Time evolution of a scalar field on one moving wall 
for V — 0.4, A = 1.0. The value of the scalar field is given by 
^w{t) = $(t, yw(i)), where J/w(i) is the position of the wall 
and r is the proper time on the wall. 
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FIG. 5: Time evolution of a scalar field "l>w(''") on one moving 
wall for V = 0.2, A = 1.0. 



21 and [S| One may wonder whether this osciUation is re- 
ahstic or not. This oscillation, however, turns out not 
to be a numerical error but a real oscillation of the do- 
main wall. In Appendix [nj using perturbation analysis 
we show there is one stable oscillation around the kink 
solution $x(y)- We expect that the oscillation is excited 
by the collision. In fact, the amplitude of the oscillation 
increases as the incident velocity v increases. In large 



velocity limit (u>0.6), we find $^ w 0.18(7 — 1), where 
$00 is the amplitude of the post-oscillation. 

Since the scalar field on the domain wall oscillates as 
^w ~ 'i'oo cos ar after the collision, our wave equation 
(|3.7f) would be rewritten as 



ipk 



2^ °° 



(1 + cos2CTr) 



^fc = 



(3.20) 



where a — •\/3/2A"'"/^ is the eigenvalue of the perturba- 
tion eigen function, it is so-called Methicu equation. 

From this equation, we may wonder whether we can 
ignore this oscillation when we evaluate particle pro- 
duction rate. In fact, when we discuss a preheating 
mechanism via a parametric resonance with a similar 
oscillating behavior|3J|, in conventional cosmology^ as 
well as brane world cosmology |2a|- It is known that 
for Methieu equation, there is an exponential instability 
V^fc c)c exp(/Xfct) within a set of resonance bands, where 
/ifc = g$^/8. This instability corresponds to an expo- 
nential growth of created particles, which is essential in 
the preheating mechanism. In order to get a successful 
particle production by this resonance instability, how- 
ever, we have to require large value of ^^$^. However, in 
the present simulation, it is rather small, e.g. $00 '^-^ 0.1 
for the case of u = 0.4. Hence, we may ignore such a par- 
ticle production by a parametric resonance in the present 
calculation. However, if the incident velocity is very fast 
such as a speed of light, we may find a large oscillation. 
Then we could have an instant preheating process by the 
domain wall collision. We also wonder whether or not 
the standard reheating mechanism due to the decay of 
oscillating scalar field is effective. In this case, we have 
to evaluate the decay rate T^ to other particles. Since 
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FIG. 6: Time evolution of a scalar field "1>w(t) on one moving 
wall for V = 0.4, A = 10. 
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FIG. 7; Time evolution of a scalar field 'l>w(''") on one moving 
wall for V = 0.2, A = 10. 
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TABLE I: The number and energy densities (n and p) of 
created particles for the typical values of the coupling g, the 
incident velocity v and the self-coupling A. d — \/2/\ and 
Nh denote the width of the wall and the number of bounces 
at the collision, respectively. 



r^ ex g"*, we expect that the reheating temperature is 
proportional to (f' , which is enough small to be ignored. 
Note that there is another factor which reduces the decay 
rate in the case that the potential is not a spontaneous 
symmetry breaking typei34!|. 

In what follows, we ignore the creation due to post- 
oscillation stage. Hence, we just follow the procedure 
shown in the previous subsection. Using the evolution 
of the scalar field $Wi we calculate the Bogolubov coef- 
ficients ttfc and /9fc . In Table U we show the results for 
three different parameters; v (the incident velocity), A 
(the self-coupling constant of the scalar field), and g (the 
coupling constant to a particle ip). 

From this Table, we find the following three features: 



(1) The produced energy density p depends very much 
on g. We study two cases with g = 0.01 and 0.1. The 
energy density for g — Q.\ is 10'' times larger than that 
for g = 0.1, which means that p is proportional to g*. 

(2) The energy density p for v — 0.2 is twice larger than 
that for V = 0.4. It may be so because the bounce occurs 
twice for v = 0.2, while once for v = 0.4. 

(3) The energy density is less sensitive to A. 

We also investigate for several different initial veloc- 
ities because the coUisional process is very sensitive to 
its incident velocity. We analyze many cases with two 
bounces, with three bounces, with four bounces, • • • in 
the range v — 0.2-0.25 as shown in Appendix B. Using 
those numerical data, we also evaluate the number and 
energy densities of the particles created at the collision. 
The results are summarized in Table Hll 

From this Table, we confirm the above three features 
(l)-(3). In particular, it becomes more clear that the 
energy density is proportional to the number of bounce 
Nb. 

We can summarize our results by the following empir- 
ical formula 
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(3.21) 
(3.22) 



If this energy of the particles is thermalized by inter- 
action and thermal equilibrium state is realized, we can 
estimate the reheating temperature by 
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(3.23) 



where ^off is the effective number of degrees of freedom 
of particles. Hence we find the reheating temperature by 
the domain wall collision as 



Tn= ^ 
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In order to see more details, in Figs. |Sland|ni we show 
a spectrum of the produced particles of number density 



n, I.e. 



n = j dk Uk with Uk = 47r|/?fc| k 
/o 



(3.25) 



The spectrum nk is well fitted as a gaussian distribution 
as 



Hk ~ inAe ^''0 



(3.26) 



where fco = 0.73 and A = 3.12 x 10"^ for Fig. 10 and 
fco = 2.04 and A = 3.43 x 10"^ for Fig. 11, although 
there is small deviation partially. These parameters can 
be described by physical quantities as fcp ~ 1/d and 
A Ri ^od^g'^. The reason is well understood, fco ~ 1/d 



means that the typical wave number is given by the width 
of the scalar field when domain walls collide (see Fig. ^. 
As for /3, it corresponds to the "reflection" coefRcient of 
the "potential" given by Figs. I4I7I and then it will be 
proportional to the coupling constant g^, and the reflec- 
tion rate (|/3p) will be related to the potential depth $o 
and the square of the width d? . 

Integrating the fitting spectrum 13.26|l . we obtain 

POO 

n = Ukdk = (27r)3/2$orf5^ ~ 25^5^ , (3.27) 

/•oo 

p = / UkUJkdk = 47r$o5^ « 20^'' , (3.28) 

Jo 

which values are exactly the same as those obtained by 
numerical integration in the case with one bounce (see 
Eqs. H3.21|l and H3.22|l ). We expect that they are en- 
hanced by the factor Nb when we find iVf, bounces at the 
collision. 

Therefore, although we obtain the particle creation nu- 
merically, the result is easily understood and summarized 
by a simple formula. 

IV. SUMMARY AND DISCUSSION 
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FIG. 8: Spectrum of the created particles at the collision. We 
plot log |/3fcpfc^ with respect to k for v — 0.4, A = 1.0, g = 
10^^. The gaussian distribution is plotted by a dotted line, 
which gives a good approximation for k < 3. In the small 
box, we enlarge the low frequency region {k < 1) to see the 
deviation from the gaussian distribution. 



We study the particle production at the collision of two 
domain walls in 5D Minkowski spacetime. This may pro- 
vide the reheating mechanism of an ekpyrotic (or cyclic) 
brane universe, in which two BPS branes collide and 



Nb 


V 


A 


d 


n 


P 


2 


0.225 


1.0 


1.414 


7.03x10-^ 


3.72x10"^^ 


10 


0.447 


2.21x10"'^ 


3.71x10"^^ 


0.238 


1.0 


1.414 


7.08x10-^ 


3.78x10"'^ 


10 


0.447 


2.23x10"'^ 


3.78x10"'^ 


3 


0.2062 


1.0 


1.414 


1.10x10^'' 


6.07x10"^ 


10 


0.447 


3.45x10"'^ 


6.06x10"'^ 


0.2049 


1.0 


1.414 


1.09x10^'' 


6.01x10"'^ 


10 


0.447 


3.43x10"^^ 


6.01x10"'^ 


0.2298 


1.0 


1.414 


1.10x10"'' 


6.04x10"'^ 


10 


0.447 


3.43x10-'^ 


6.02x10"'^ 


0.22933 


1.0 


1.414 


1.09x10"'' 


6.03x10"'^ 


10 


0.447 


3.44x10"^^ 


6.01 xlO"'' 


4 


0.229283 


1.0 


1.414 


1.47x10"" 


8.10x10"^^ 


10 


0.447 


4.61x10"'^ 


8.09x10"'^ 


0.2292928 


1.0 


1.414 


1.47x10"" 


8.16x10"^^ 


10 


0.447 


4.62x10"'^ 


8.17x10"'^ 



TABLE II; The number and energy densities (n and p) of 
created particles with respect to the number of bounces A^;,. 
V and d are the incident velocity and the width of the wall, 
respectively. We set g = 0.01. 



evolve into a hot big bang universe. We evaluate a pro- 
duction rate of particles confined to the domain wall. 
The energy density of created particles is approximate 
as p « Air^og'^Ni, where $o is the maximum amplitude 
of $-W) ^b is the number of bounce at the collision and 
g is a coupling constant of a particle to the scalar field 
of a domain wall. If this energy is converted into stan- 
dard matter fields, wc find the reheating temperature as 
Tr w 0.88 X g N^^'^ig^s/lOOy^^^. We find that the 
particle creation is affected more greatly by a coupling 
constant g than the other two parameters v and A. The 
initial velocity changes the collision process, that is the 
number of bounce at the collision, but this is less sensitive 
to the temperature. The thickness of a domain wall d (or 
a self-coupling constant A) changes the width of poten- 
tial ^1,{t) of a particle field (ip), and it changes a typical 
energy scale of created particles, which is estimated as 
w - 1/d. 

In order to find a successful reheating, a reheating tem- 
perature must be higher than 10^ GeV, because we wish 
to explain the baryon number generation at the electro- 
weak energy scale 2] . Since Eq. H3.24|l is written in the 
following form; 
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(4.1) 



we find a constraint on a fundamental energy scale ttt,^ as 



m, 



,>1.1 X 10^ GeV for g = 10^^ 



and TO^>1.I X 10 
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FIG. 9: Spectrum of the created particles at the coUision. 
We plot log |/3fc pfc^ with respect to A; for u = 0.4, X= 10,g = 
10~^. The gaussian distribution is plotted by a dotted line, 
which gives a good approximation for fc < 2.5. In the small 
box, we enlarge the low frequency region (fc < 1) to see the 
deviation from the gaussian distribution. 



for g = 10 ^, which are slightly larger than TeV scale. 
Here we assume gcs = 100 

In the present work, we consider the 3-dimensional do- 
main walls in the 5-diniensional Minkowski space and 
show that the particle production at the two-wall colli- 
sion may provide a successful mechanism for reheating in 
the ekpyrotic universe. In a string/M theory, however, 
we expect higher dimensions, e.g. 10 or 11. If we com- 
pactify it to the effective 5-dimensional spacetime, our 
work can be applicable to such a mode. Moreover, if we 
discuss a collision of a p-dimensional walls (branes) in 
{p + 2)-dimensional spacetime, our approach can also be 
extended. 

In this paper, we have not taken account of effects from 
a background spacetime. We are planning to study how 
such a generalization affects the present results about 
particle creation at the collision. 
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APPENDIX A: NUMERICAL METHOD 

For our numerical analysis of the domain wall colli- 
sion, we solve the partial differential equation H2.3|) on 
discrete spatial grids with a periodic boundary condi- 
tion. The scalar field on the grid points is defined by 
$„(i) = $(y„, t), where y„ = nAy, for n = 1, 2, . . . , iV. 
We use the fourth-order center difference scheme to ap- 
proximate the second spatial derivative [3g as 



52$, 



1 



dy^ -12(A,)^[-*"- + l^'^"--^«'^" 
+ 16$„+i - $„+2] + 0((Ay)4) 



(Al) 



This leads to a set of N coupled second-order ordinary 
differential equation (ODE's) for the $„, i.e. 






1 



12(Ay)2 ' 
+16$„+i 



- $„+2 



16$„_i - 30$, 
-A$„($2_i) 



(A2) 



The ODE's IIA2p are solved using a forth-order Runge- 
Kutta scheme, and so our numerical algorithm is accurate 
to fourth order both in time and in space, with error of 
0{{Ay)'^) and 0{{AtY). For the boundaries, we set the 
left and right grid boundaries at y^ = —40 and yj^ — 
+40, and impose the condition as <&(?; = VL^i) = ^{v — 
yR,t) = —1 . The grid number is TV = 8000 with the grid 
size of Ay = 1.0 x 10^^. The initial position of a wall j/o 
is set by yo = Ivl + {vr - Vl + 1)/3| = 13, equivalently, 
one-third of the numerical range. For time steps, we set 
At = 0.7 X Ay. 

As for the particle production process, we have to solve 
the second-order ordinary differential equations H3.7|) for 
each wave number k. By using the fourth-order Runge- 
Kutta scheme, we solve them for the wave number k of 
< A: < 100 with the width Afc = 1.0 x lO'^. We esti- 
mate |/?fcP in the equation (|3.13|) . Defining the functions 
like these 



Wi 


= i'ki^l - ^k'^k , 


(A3) 


W2 


= ^ki^l + ^liPk , 


(A4) 


W3 


= Ark , 


(A5) 


we use the formula 






15 P ^^'' 


Wi/2iLu)^ + {W2/2lu)^ 


(A6) 


\Pk\ 


AW3 


to evaluate \Pk\- 







APPENDIX B: NUMERICAL EXAMPLES OF 
SEVERAL BOUNCES AT THE COLLISION 

We depict some numerical examples which show sev- 
eral bounces at the collision. First we show two typical 
examples in Fig. ^| The figures show the behaviors 




(a) 



(b) 



FIG. 10: Time evolution of <l> field at y = are depicted for 
(a) w = 0.4 A = 1.0 and (b) v = 0.2, A = 1.0. The bounce 
occurs once for a large velocity (a), while Two bounces are 
found for the slower velocity. 



a domain wall structure which is excited by the collision. 

From Fig. ^| we find there is one bounce for v = 0.4, 
while a bounce occurs twice for v — 0.2. In fact, the 
results are v ery much sensitive to the incident velocities 
as shown in [29J. 

Here we present several examples to show how the be- 
haviors of the scalar field depend on v, which confirm the 
previous works. 

In Figs. EKa-) and (b), we show the case for v = 0.2062 
and V = 0.2049, respectively. We find three bounces 
at the collision. Four-bounce solutions are depicted in 
Figs. CKa) (v = 0.229283) and (b) {v = 0.2292928). 
These calculation shows that the detail collisional process 
is very sensitive to the incident velocity. 



APPENDIX C: PERTURBATIONS OF A 
DOMAIN WALL 
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(a) 



(b) 



FIG. 11: Time evolution of "I> field at y = are depicted for 
(a) V = 0.2062 A = 1.0 and (b) v = 0.2049, A = 1.0. Three 
bounces are found. 



of the scalar field $ at y = with respect to t. Ini- 
tially, when two domain walls locate at large distance, 
the value of the scalar field at y = is 1. Then two 
walls approach and collide. At this point the value of 
1$ — 1| increases. After the collision, it again decreases 
to an initial value. We find some small oscillation around 
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(a) 



(b) 



FIG. 12: Time evolution of <E> field at y = are depicted for 
(a) V = 0.229283 A = 1.0 and (b) v = 0.2292928, A = 1.0. 
Four bounces are found. 



We show that the oscillation we find in Sec IIII Bl is a 
proper oscillation around a static stable domain wall. To 
show it, we perturb the static domain wall solution H2.4|l 
as 



$ = <^K{y)+S<^{t,y) 



(CI) 



In this appendix, we use the dimensionless variables 
rescaled by rj. Substituting this into Eq. H2.3|l and lin- 
earizing it, we obtain 



(5$ - (5$" + A(3$^ - 1)5$ = 



Setting 



(5$ 



'F{y) , 



(C2) 



(C3) 



and introducing new variable y = tanh(y/(i), we rewrite 
Eq. (EH) as 



d^F 



dF 



(l-r)-Trj-2y— 
dy dy 



3- 



2-(tVA 



1 



y' 



F = . (C4) 



The solution is given by the associated Legendre function. 
Imposing the boundary condition (F —^ as y —^ cxd), we 
have two regular solutions; F{y) — P-^iy) with cr = 
and P^{y) with a = ^/zj2\^l'^ . 

For the former mode, it just corresponds to a boost 
of a kink solution in the y direction, because we find 
F(y) K, ^Kiy) — ^K{y — dy). The latter one is the os- 
cillation mode which we find in Sec 111 Bl In fact, taking 
the average over 10 cycles in the oscillation after the col- 
lision in Fig. ^ we obtain the mean angular frequency 
is about 1.17, which is very close to cr = ■\/3/2 r^ 1.22. 
The ratio is about 0.96. We also evaluate the angular 
frequency for other cases. We find 1.33 = 1.09(T for Fig. 
3.67 = 0.95CT for Fig. and 4.18 = 1.08cr for Fig. 
From the figures, we also find that the amplitude of oscil- 
lation gets larger as the incident velocity is faster. This 
is because the excitation energy of a wall at the collision 
will be large for a large velocity. 
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We then conclude that the osciUations after the colh- 
sion of two domain walls is the proper oscillations around 



a stable domain wall. 
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